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Abstract 

The main objective of this paper is the construction of the solution of 
an impulsive stochastic differential equation, subject to control conditions 
in the pulse times, and to give sufficient conditions for the pulse times to 
be random variables with finite expectations. Such equations are useful 
in modeling diverse phenomena, such as biological control and pressure 
regulating mechanisms. The article ends with an application to fisheries. 


1 Introduction 


Stochastic differential equations, SDEs, have been an important tool in the de¬ 
scription of many phenomena in different areas of knowledge. The main results 
as to their existence, uniqueness, and qualitative properties are found in the 
classical literature, such as 1|8|6 and in recent years, a line of research has 


been exploring problems related to stochastic differential equations with pulses 
(see |5|10|11| ). 


This paper focuses on determining the finiteness of the expectation of the 
pulse times that occur when the solution process of a SDE reaches a control 
function s{t), producing a pulse that sends the process to a second control func¬ 
tion q{t). 


To this end, consider a one-dimensional Brownian motion Bt defined over 
a stochastic basis (fl, T, (T)t, P), two real and continuous functions q and s, and 
a stochastic process X{uj,t) : (fl, R+) —>■ R solution of a SDE of the following 
type: 

dA(f) = f{X{t))dt -I- aX{t)dBt, 

with Ai(0) = g(0) and r = inf{t > 0;A(t) = s{t)}. In this time, the process 
will be sent to the value X{t^) = qij). 

The functions q and s represent the curves of the lower and upper control 
functions, respectively. 


2 The Model 

Let q and s be two continuous and positive functions such that q{t) < s{t) for 
all t > 0, let / S C^(R), and let Bt be a standard Brownian motion. Then the 


1 




model can be written as 


dX{t) 


f{X{t))dt + aX{t)dBt 

q{Tk) 

inf|t>rfc_i; X(t) = s(t)| 




To = 0 


( 1 ) 



This figure shows a possible path of evolution of the system proposed. The 
times ti,T 2 and T 3 correspond to the first pulse times of the system. 

The main goal is to prove that it is possible to obtain a behaviour as in Figure 
1, that a local solution of the equation dX{t) = f{X{t))dt + aX{t)dBt exists 
on the stochastic interval ]rfc_i,Tfe], and that the pulse times Tk are random 
variables with finite expectations. 


3 HYPOTHESES 

(A) Suppose that / e C^(IR+) and that there exist positive constants a,/3 and 
S such that 0 < a < (3 and 

Oix < f{x) < fix, X € [0, S]. 


(B) 


Assume that the positive constants a and a satisfy a > 


a 


1 


2 ■ 


(C) Assume that q and s are two continuous functions such that 0 < q{t) < 
s{t) < S for all t>0. 
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NOTE: Hypothesis (A) implies that 



• /(O) 

then 


n -f 
0 , if — 

ax 


exists, 


21 = 0 


• a< /'(O) < p. 


Figure 2 shows the relation of Hypothesis (A). 


Proposition 1. Suppose that f satisfies Hypothesis (A) and that xg g]0, S'[ 
and a are positive constants. Then there exists a unique stochastic process X (t) 
such that 


( dX{t) = f{X{t))dt + aX{t)dB{t) 

\ A(0) = Xg 

on the stochastic interval [0, r] where r = inf {t > 0; 0 < A(t) < S}. 


( 2 ) 


Proof. Hypothesis (A) implies that / is a Lipschitz and sublinear function on the 
stochastic interval [0,r], where r = inf {t > 0; 0 < X{t) < S}. Therefore the 

application of the classical theorem of the existence of a solution to a stochastic 
differential equation with local Lipschitz coefficients (see BaldP or 0ksendaPl) 
implies that ([^ has a unique solution process on the interval [0,t]. □ 

NOTE: Under Hypothesis (A), the process X{t) is a stopping time with con¬ 
tinuous paths, therefore 


r = inf {t > 0; 0 < X{t) < S} = inf {t > 0; 0 < X{t) = S} . 

Proposition 2. Suppose that f and a satisfy Hypothesis (A). Then the stochas¬ 
tic process X (t) that is the solution of (§ satisfies 


< X{t) < Xi^{t), (3) 

where A„(t) = and Xp{f) = 

Proof. Note that if f{x) < g{x) with x G I, then xi{t) < X2{t), where xi 

and X2 are the solutions of the equations x' = f{x) and x' = g{x) with 

a:i( 0 ) = 2:2(0) = Xg. 


The fact that ax < f{x) < Px, x e [0, S] implies that the processes Xa{t), 
X{t) and Xp{f) are solutions to the equations 
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(i) dX(t) = aX{t)dt + aX{t)dB{t), 

(a) dX(t) = f{X{t))dt + axlt)dB(t) and 
(in) dX{t) = (]X{t)dt + aX{t)dB{t) 

respectively, with Xa(0) = X(0) = -^,s(0) = xq £]0, S'[ 


satisfying Xa{t) < X{t) < Xp{t) on the stochastic interval t = inf {t > 0; 0 < X{t) < S}. 


Furthermore, the explicit solutions of (i) and (in) correspond to the geometric 

Brownian motions X^ = a;oe(“ ^ 2 

respectively (see 0ksendal,l^ p. 63). So the result is obtained. 


□ 


Note: Proposition implies the positivity of the solutions. If Hypothesis (B) 
is not satisfied, and in particular if a < —, then 


0 < X{t) < Xi 3 {t) —>■ 0 where t —)■ 00 a.s. 


Theorem 1. If f and a satisfy Hypotheses (A) and (B), then the stopping time 
T = inf {t > 0; 0 < X(t) = S'} is a random variable with finite expectation. 

Proof. From Propositionwe have Xa{t) < X(t), Vt > 0, therefore t < Ta, 
where r = inf {t > 0; 0 < X(t) = Sj and = inf {f > 0; 0 < Xa(t) = S}. 


From the equation Xa(t) = S we have that 

.2x r s 


- A- r + — 


Xo 


and therefore 


E(r„) = 


ln(^ 

xq 


( 4 ) 


By hypothesis, xq < S and a > —, and so 0 < E(tq,) < 00 , completing the 
theorem. 


□ 


NOTE: Similar arguments can be used to deduce that E(t/ 3 ) < E(r), where 
Tp = inf{t >0; 0 < Xp(t) = S}, therefore the expectation for r is bounded 

by 


E(r/3) 



< E(r) < 



E(t„). 


The following figure shows the relation between r^, r and rp. 


(5) 
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A study of the first exit time of diffusion, can be seen in [^, but with 
assumptions and different approach to this work 


4 The pulses and the control 


The interest of this paper is to bring the process from a state on the curve q{t) 
to a state on the curve s{t) and return it, instantaneously, to a state in q{t), as 
is shown in Figure 1. Therefore consider the equations before the first pulse. 


J dA(f) = f{X{t))dt + (jX{t)dBt 
\ A(0+) = q(0). 


and the random variable Ti 


= inf > 0; 


X{t) = 


( 6 ) 


Corollary 1. Under Hypotheses (A), (B) and (C), Equation Q has a unique 
solution on 0 < X(t) < S, and ti is a random variable with finite expectation. 


NOTE: Due to the continuity of X{t) and s{t), the random time ri is well 
defined. 


Proof. Due to Hypothesis (C), A'(O) > 0 and A'(O) < s{t) < S, Vf, it is possible 
to use Proposition]^ Therefore, there exists a unique stochastic process Xi{t) 
solution of Equation ([^ on the interval ]0,ri]. The fact that sit) < S implies 
that Ti < r, Vw S D, where r is the stopping time defined in Proposition 
taking xq = (?(0). Therefore E[ti] < E[t] < oo. □ 

Corollary 2. Let T > 0. The equation temporally displaced by T, 

I dX(t) = f(X(t))dt + aX(t)dBt 
{ X(T+) = q{T), 
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has a unique solution on 0 < X(t) < S, and the random variable defined by 
T 2 = inf > T; X(t) = s(t)'^ has a finite expectation. 

Proof. The change of variables u = t — T transforms Equation Q into 

f dX(M) = f{X{u))du + aX{u)dBi^^ /o', 

1 X(0+) = q{T), ^ ’ 

where B^u'^ = B^ — Bt is a new standard Brownian motion and r 2 = T + r, 
with T the stopping time defined in Proposition taking xq = q{T). 

Therefore, it is possible to apply Corollaryj^with initial condition X(0+) = q(T) 
and Brownian motion B^^'>, thus obtaining a stochastic process X{u) solution of 
Equations Q with u s]0, rfi that is equivalent to the solution X(t) of Equations 
Q with t S]T, r 2 ]. Moreover E[r 2 ] = T + E[r] < oo. □ 

Note that between E[t 2 ] and E[t] there is no relationship of order. 


Theorem 2 (Main results). Under Hypotheses (A), (B) and (C), there exists 
a stochastic process X(t) that is the unique solution, in the a.s. sense, of the 
system Q, and the pulse times t/- are random variables with finite expectations. 

Proof. Using Corollaryand Corollaryit is possible to obtain Xift), a so¬ 
lution of Equation § on the interval ]0,ri], and X 2 {t), a solution of Equation 
0 with initial condition = q{Ti) on the interval ]Ti,r 2 ]. 

Then the process Xlt) = Ai(t)l]oi-i](i) is a solution of the system Q on the 
stochastic interval ]0ti] and X{t) := X2(t)ll]i-i,r2](0 is a solution of the sys¬ 
tem Q on the stochastic interval ]ti,T2]. Moreover, X{t) = Ai(t)l]o,-p (t) -|- 
A 2 (t)l]T-i,T 2 ](^) is a solution on the interval ]0,T2]. Therefore, using the same 
argument recursively, we find 


k^l 

a global solution of the system Q. Its uniqueness a.s. is due to the uniqueness 
on each interval. 

The expectation of the stopping time = inf > Tk-i\ X(t) = s(t)| is 
given by the following recurrence equation. 

E[Tfe] = E[Tfc_i] -I- E[ATfe] (10) 

where Arfc is the stopping time defined in Propositionj^taking Xq = q(Tk-i). □ 

The random variable At^ = — Tk-i corresponds to the time between two 

successive pulses, called the fcth timeout. 
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5 Asymptotic behaviour 

In order to have more precise information about the asymptotic behaviour of 
the process, it is necessary to impose additional conditions on the functions q 
and s. More precisely, if q is increasing and s is decreasing, or vice versa, then 
it is possible to improve the estimates for the timeout Ar^. 

Proposition 3. Assume that Hypotheses (A), (B) and (C) are satisfied and 
further, let s be a decreasing function and q be an increasing function. Then 
E[Arfc] is a decreasing sequence. 

Proof. Consider the auxiliary times Tq and defined by 

Tq = inf > Tk-i; Xk{t) = g(Tfc)| y = inf > r,; Xk{t) = s(Tfc+i)|. 

( 11 ) 

Note: (Xk, t^) is the solution to the system Q, and the fact that q is increasing 
and s decreasing implies Tf^-i < Tq < Tg < tj.. 

The following graph describes the situation: 



Arfe = {Tq - Tk-l) + {Ts - Tq) + (rfc - Ts) 


( 12 ) 


However, the random variable 


corresponds to the timeout of Equation 


0 with initial time T = Tq, initial value q{Tk), and final value s(rfe_|_i), and is 
therefore a random variable distributed identically to the timeout Ar^+i. Now 
taking expectations in Equation (121, we have 


E[Arfc] = E[t, - Tk-l] + E[ts - Tq] + E[Tk - Ts] 

= E[t, - Tfc_l] + E[ATfe+l] + E[Tfe - Tsj. 

Therefore E[ATfc] — E[ATfe_|_i] = E[Tq — Tk-i] + E[Tfe — Ts] > 0 , completing the 

proof. □ 
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Proposition 4. Assume that Hypotheses (A), (B) and (C) are satisfied, and 
that s is an increasing function and q is a decreasing function. Then E[Arfc] is 
an increasing sequence. 


Proof. The proof is similar to that of the previous proposition. In this case, 
the temporal order is Tk-i < Tk < Tq < Tg < Tk+i^ with Tq and Tg defined 
analogously to Equation The k + 1th timeout is decomposed as follows. 


At/c+i = {Tq - Tfc) + (Ts - Tq) + (t^+i - Tg) 


(13) 


Now, the random variable Tg — Tq corresponds to the timeout of Equation 0 
with initial time T = Tq, initial value q{Tk-i), and final value s(rfc). It is 
therefore a random variable identically distributed to the timeout Arfc. Now 
taking expectations in Equation (121, we have 


E[Arfc+l] = E[t, - Tfc] + E[Ts - Tq] + E[Tk+l - Tg] 

= - Tfc] + E[Arfc] + E[Tfe+i - Ts]. 

Therefore, E[Arfc+i] — E[ATfc] = E[Tq — Tfc] + E[Tfc+i — Tg] > 0, which completes 
the proof. □ 


Note that if q and s are monotone, then under Hypothesis C, lim q{t) = 

>-oo 

Q S]0, S'] and lim s{t) = S g] 0, S[ with 0 < Q < S < S. This allows the 

t—>-oo 

following theorem. 

Theorem 3. If Hypotheses A, B and C are satisfied, if q is increasing and s 
is decreasing, or if q is decreasing and s is increasing, then the expectation of 
the timeout Arfc converges on the interval [a/ 3 ,aa], with 


Proof. 



and Qa 



1. If s is decreased and q is increased then, by Proposition]^ the sequence 
of the expectations of a timeout, E[ATfc], is a decreasing sequence. In 
the expression for the expectation of t^ in Equation (j^, the initial and 
final values can be replaced and updated at time fc by Xg —> q(Tfc_i) 
and S —> q{Tk-i). Therefore the expectation of the A:th timeout Arfc is 
bounded by 


1 



In 


/ s(Tfc_l \ 

\q{Tk-i)) 


and taking limits, 


< E[ATfc] < 


1 



In 


/ s(Tfc-i \ 

\qiTk-i)J ’ 


(14) 
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< lim E[ATfc] < 

fc—>-oo 



(15) 


and so E[ATfc] is a decreasing and bounded sequence, therefore it has a 
limit in the interval [ap, Oa]. 











2. The case where s is increasing and q is decreasing is analogous to the 
previous case, but then E[Arfc] is an increasing sequence. 


□ 


6 Applications 

1. Fishery with fixed quota: We consider a fishery resource with long 
periods of closure. When the amount of resources reaches a certain value 
S, fishing is allowed for a very short period of time with a total fishing 
quota C, and then the closure of the resources begins again. The evolution 
of the resources follows a logistic growth law in a random environment. 
This can be modeled by the following impulsive system. 


dX{t) = r (^1 - X{t)dt + aX{t)dBt t e]Tk,Tk+i] 

< X{r+) = S-C 

Tk = inf|t>rfc_i; A(t) = S'!, tq = 0, 

(16) 


where r is the intrinsic growth rate and K the capacity of the environ¬ 
ment. It must be verified that Hypotheses (A), (B) and (C) are met. 


The population dynamics in a random environment has been extensively 
studied by many authors; a complete compendium of the main results can 
be found in [^. Applications to fisheries have been studied in [^. The 
regulation of fishing and fishing quotas in a random environment have 
been addressed in [^. Impulsive models of a fishery have been treated in 
and 

• Hypothesis (A): In this case, f{x) = r ^1 — x, the constant S 
may be chosen as any positive value less that K. Thus, the problem 
makes sense. 
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fiS) r 

The constant = f'{0) = r and the constant a = —^ = —{K—S). 

o J\ 

Hypothesis (B): The constant a can be anything such that < 2a. 
Hypothesis (C): In this case, s{t) = S and q{t) = S — C. 


Due to the fact that q and s are constants, the timeouts are indepen¬ 
dent and identically distributed random variables, therefore E[Afc] = E is 
constant on the interval [ 0 ^, 0 ^], with 


fl/3 — 2 

r — 


In 


S 

S-C 




s 

S-C 
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2. Fishery with programming of total closure: The context of this 
application is a resource that will have a total prohibition on its being 
fished, but during a period there will be allowed a quota, progressively 
smaller, only for short periods of fishing. As in the previous case, the 
fishing is permitted only when the resources reach the value S, but the 
quota is a fraction of the amount of the previous fishing, depending on 
the timeout: this is 


C{n) 

C(rfe_i) 


where 7 £]0,1[ is the fraction of fishing and T is a temporal scale factor. 
Note that if is constant, then C(rfc) is a geometric progression. A 
reasonable choice for C is C'(rfc) = , as in the previous example, 

q{t) = S — C, and therefore q{t) = S (^l — 7 ^+T^ and the model can be 
written as 


dX(t) = r(^l-^^^^X{t)dt + crX{t)dBt t G]Tk,Tk+i] 

< A(r+) = 5 (1-71+^) 

Tk = inf|t>rfc_i; X(t) = S'!, tq = 0. 

(17) 


As 7 e]0,1[, then g(0) = S'(l — 7 ) S ]0, S[ and 


As a consequence of Equation (151, one obtains 


lim q{t) = S. 

k—¥oo 

that lim E[ATfc] = 0 . 

A;—>-oo 
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